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ABSTRACT 

The rate of turbulent spherical accretion onto a compact object might be much 
smaller than the Bondi rate. It is suggested that the rate of accretion onto Sgr A* is 



Q H . much smaller than the Bondi rate. 



Subject headings: accretion — turbulence 



I/"") ' 1. The rate of spherical accretion 

The rate of adiabatic inviscid spherical accretion onto an object of mass M is given by the 
Bondi formula 

OO . '2 

M Bon di ~ PC S R A , (1) 

i-Q ■ where p is the density of the accreting g cLSj Cg IS the speed of sound, and Ra is the accretion radius 

Q-f 

6 1 R A ~ GMc~ 2 . (2) 

We show that turbulent heat conduction reduces the rate of accretion onto a compact object 

x. *-*-(«![)"• (3) 

where Rs -C Ra is the radius of the compact object, and a is a model-dependent positive parameter. 
This formula is derived in the next section. Turbulent viscosities are not considered in this work, but 
they should also slow down the accretion. Turbulent heat conduction hinders spherical accretion 
because it saves part of the heat from being accreted onto the object. The accretion flow becomes 
hotter, and hot gas accretes slower. 

To illustrate the astrophysical applications of this result consider the case of Sgr A*, which is 
believed to be a 2.5 x 1O 6 M0 black hole located at the center of our Galaxy; its accretion luminosity 
might be L ~ 10 37 erg/s (Genzel et al 1994). When stellar winds in the vicinity of the galactic center 
collide and shock, they produce a gas of c s ~ lOOOkm/s and p ~ 10- 20 g/cm 3 (Coker & Melia 1997). 
The Bondi accretion rate of this gas is Meondi ~ 10 21 g/s. Then MBon&iC 2 ~ 10 42 erg/s is five orders 
of magnitude higher than the actual luminosity. We (Gruzinov 1998, Quataert &; Gruzinov 1998) 
and others (Meszaros 1975, Blandford 1998) have argued that such low radiative efficiencies are 
unrealistic. It is then natural to assume that Sgr A* accretes at a much smaller rate (Blandford & 
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Begelman 1998). If we use the estimate (3), and arbitrarily assume a = 0.4, the estimated accretion 
rate is reduced by a factor of 100. The radiative efficiency is then 0.1%. This radiative efficiency is 
still small, but it might be reasonable because electrons might be heated much less than ions, and 
even hot electrons might radiate inefficiently when the plasma is rarefied. 



2. A solvable model 

Here we derive equation (3). We assume a spherically symmetrical gas inflow. The flow is 
taken to be inviscid but non-adiabatic, i.e. we include turbulent thermal conduction and neglect 
turbulent viscosities. This artificial assumption is made to simplify the model. No real calculation 
is possible anyway, since we do not understand the nature of MHD turbulence in the accretion flow. 
Our sole purpose is to demonstrate that the rate of turbulent spherical accretion might be much 
smaller than the Bondi rate. 

The gas density p, the inflow (positive) velocity v, and temperature T depend on the radial 
coordinate r only, and satisfy the stationary equations of continuity, Euler's, and the thermal 
conduction 

pvr 2 = J, (4) 
vv' = -p- 1 (pT)' -r~ 2 , (5) 

PTV {~p ~ It) = '~V 2 T7- (6) 

Here J = M/Att, the prime denotes the radial derivative, GM = 1, m p = 1. We assume that the 
thermal conductivity k is given by the "Shakura-Sunayev" formula 

k = aprv, (7) 

where a ~ 1 is a dimensionless positive constant. 

The boundary condition at the surface of the compact object is unknown, but it is irrelevant. 
For concreteness, assume T + (3rT' = at r = Rs, where (3 is a dimensionless constant. At r = oo, 
we assume T = 1, so that Ra ~ 1- The boundary condition for p is also irrelevant, because the 
system of equations is invariant under p — > Xp, J — > A J. 

Our aim is to find a maximal possible value of m = J/p(oo), for which a smooth solution of 
the system exists. The ratio M/Meondi is equal to this maximal value. Obviously m ~ 1 if a ~ 1 
and Rs ~ 1. We need to find out how m depends on Rs when Rs <C 1. 

Scale out J by p — > Jp, plug (4) and (7) into (6): 

pvr 2 = 1, (8) 



vv > = -p-\pTy - r - 2 , 



(9) 
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T ( — — ^777") = a(rT')'. (10) 



p 2 T 

Using (10), integrate (9) and obtain a system of two first-order differential equations 

115 ,15 

2VV + 2 T + arT =~r + r (U> 

p V r P / r r°p z 
We need to find how the maximal possible value of m = l/p(oo) (for which a smooth solution of 
the system exists) depends on R$ <C 1. For small r, we can neglect 5/2 in the right hand side of 
(11). Then, denoting r = rT, f = (r 3 p 2 ) -1 , x = logr, we obtain 

cLt 

2a— = 2 - (5 - 2o)t - f, (13) 

df_ = f 2 + 2dr/dx-5T-f 
dx t — f 

This system has a stable equilibrium point r = 2(5 — 2a) -1 , / = 0. In the vicinity of this point, 

= -2a/. (15) 

ax 

The system (13), (14) is an accurate approximation only for x < 0, that is for an x-duration 5x = 
log(l/i?s). The quantity / decreases according to (15) for an x-duration dx = log(l/i?s) — const. 
Therefore the final density is proportional to (Rs)~~ a - Since m = l/p(oo), the scaling law (3) is 
proven. To check the answer we integrated (11), (12) numerically. 
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